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Abstract 

Higher order scrambled digital nets are randomized quasi-Monte 
Carlo rules which have recently been introduced in [J. Dick, Ann. 
Statist., 39 (2011), 1372-1398] and shown to achieve the optimal rate 
of the root square mean error for numerical integration of smooth 
functions defined on [0, l] s . The key ingredient there is an interlacing 
function applied to the components of a randomly scrambled digital 
net whose number of components is a multiple of the dimension. In 
this paper, we replace the randomly scrambled digital nets by ran- 
domly scrambled polynomial lattice point sets, which allows us to 
obtain a better dependence on the dimension while still achieving the 
optimal rate of convergence. We consider weighted function spaces 
with general weights, whose elements have square integrable partial 
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mixed derivatives of order up to a > 1. We prove that the component- 
by-component construction of polynomial lattice rules can be used to 
obtain explicit constructions of good polynomial lattice rules. By first 
constructing classical polynomial lattice rules to which we then apply 
the interlacing scheme of order d we obtain a construction cost of the 
algorithm of order 0(dsmb m ) operations using 0(b m ) memory. 

1 Introduction 

In this paper we study the approximation of multivariate integrals of smooth 
functions denned over the s-dimensional unit cube [0, l] s , 

1(f) = f f(x)dx, 

by averaging function values evaluated at n points Xo, . . . , x n -i with equal 
weights, 

n— 1 

While Monte Carlo methods choose the point set randomly, quasi-Monte 
Carlo (QMC) methods aim at choosing the quadrature points in a determin- 
istic manner such that they are distributed as uniformly as possible. The 
Koksma-Hlawka inequality guarantees that such well-distributed point sets 
yield a small integration error bound, typically of order n~ 1+s (5 > 0), for 
any function which has bounded variation on [0, l] s in the sense of Hardy 
and Krause, see for instance [T8J, Chapter 2, Section 5]. Digital construc- 
tions have been recognized as a powerful means of generating QMC point 
sets [T3J 122]. These include well-known constructions for digital sequences 
by SoboP [32], Faure [15], Niederreiter [21], Niederreiter and Xing [23] as 
well as others, see [131 Chapter 8] for more information. Polynomial lattice 
point sets, first proposed in [23], are a special construction for digital nets 
and have been studied in many papers, see for example [TQ1 dU [13 [JjJl |2"U] . 
Polynomial lattice rules are QMC rules using a polynomial lattice point set 
as quadrature points. The major advantage of polynomial lattice rules lies in 
its flexibility, that is, we can design a suitable rule for the problem at hand. 

In this paper we study randomized QMC rules, that is, the determin- 
istic quadrature points are randomized such that their essential structure 
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is retained. Owen's scrambling algorithm can be used to randomize dig- 
ital nets and sequences while maintaining their equidistribution properties 
[271 128| [29] . This not only yields a simple error estimation but also achieves a 
convergence of the root mean square error (RMSE) of order n~ 3 / 2+s , for func- 
tions of bounded generalized variation. Since the estimator is unbiased, this 
can also be stated in another way, namely that the variance of the estimator 
decays at a rate of n~ 3+s . It is shown in [2] that the variance of the estima- 
tor based on a scrambled polynomial lattice rule constructed component-by- 
component (CBC) decays at a rate of n -( 2a + 1 )+ s j for functions which have 
bounded generalized variation of order a for some < a < 1. 

Here we consider higher smoothness, namely a > 1 for which we can im- 
prove the rate of convergence of the variance of the integration error further. 
The initial ideas for our approach stems from the papers [HI [7J |8] . Therein 
higher order digital constructions of deterministic point sets and sequences 
were introduced whose corresponding QMC rules achieve an integration error 
of order n~ a+s for functions with square integrable partial mixed derivatives 
of order a > 1 in each variable. An explicit construction of suitable point 
sets and sequences is the following interlacing algorithm. Let d > 1 be an 
integer and z G [0, l) ds where Zj = Zj^/b + Zj^/b 2 + • • • for 1 < j < ds. Then 
let a point x = (x±, . . . , x s ) G [0, l) s be given by 

oo d 
a=l r=l 

for 1 < j < s. Thus, every d components of z are interlaced to produce 
one component of x. To obtain a higher order digital net or sequence, one 
applies the interlacing algorithm to one of the above-mentioned digital con- 
structions. Furthermore, as shown in [5], Owen's scrambling algorithm can 
be used to achieve a convergence of the variance of the estimator of order 
n -(2mm(a,d)+i)+<5 f or a > p or d > a, this decay rate is the best possi- 
ble. For the algorithm in [9] it is important to note that one first applies 
Owen's scrambling to a point z of the digital net (or sequence) in dimen- 
sion sd and then interlaces the resulting point according to flTJ to obtain x. 
We call this method Owen's scrambling of order d, or order-ti scrambling for 
short here. In the proof of the convergence rate, it was assumed in [9] that 
the underlying point set is explicitly given by some digital (t, m, rfs)-net or 
(t, (is) -sequence. The t-value of digital (t, (is) -sequences, however, grows at 
least linearly with s, and consequently, it becomes hard to obtain a bound 
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of the variance independent of the dimension. 

In this paper, we study order-rf scrambled polynomial lattice point sets 
for numerical integration. Our strategy is to construct classical polynomial 
lattice rules in dimension ds using a suitable criterion, then apply Owen's 
scrambling to the quadrature points of the polynomial lattice rule and fi- 
nally to apply the interlacing algorithm of order d to obtain a randomized 
quadrature rule for the domain [0, l] s . We refer to such quadrature rules by 
interlaced scrambled polynomial lattice rules. The advantage compared to 
the results in [9] is the weaker dependence on the dimension and the possi- 
bility to construct the rules for a given set of weights when the integrand has 
finite weighted bounded variation, see Subsection 13.11 As in [31] , the weights 
model the dependence of the integrand on certain projections. With our ap- 
proach, we are able to obtain tractability results under certain conditions on 
the weights. 

In the next section we describe the necessary background and notation, 
namely Walsh functions, polynomial lattice rules, Owen's scrambling, and 
higher order digital constructions. In Section [3] we give a bound on the vari- 
ance of the estimator in the weighted function space with general weights 
where a function has square integrable partial mixed derivatives of order 
a > 1 in each variable. Using this bound we can introduce a quality criterion 
for the construction of interlaced scrambled polynomial lattice rules. In Sec- 
tion H] we prove that interlaced scrambled polynomial lattice rules constructed 
using the CBC algorithm can achieve a convergence of the variance of the es- 
timator of order n -( 2ram ( a > d )+ 1 )+ s . Thereafter we assume product weights for 
simplicity of exposition and describe the fast CBC construction by using the 
fast Fourier transform as introduced in [25l [26] . We show that the interlaced 
scrambled polynomial lattice rules can be constructed in order 0(dsmb m ) 
operations using order 0(b m ) memory. This is a significant reduction in 
the construction cost to previously obtained component-by-component algo- 
rithms for higher order polynomial lattice rules [I]. We conclude this paper 
with numerical experiments in Section [51 

2 Background and notation 

In this section, as necessary tools for our study, we recall the definition 
of Walsh functions, introduce polynomial lattice rules, Owen's scrambling 
algorithm, and higher order digital net constructions. 



4 



In the following, N denotes the set of natural numbers, while No de- 
notes the set of non-negative integers. The operator © denotes the digit- 
wise addition modulo b, that is, for x, y G [0, 1) with base b representations 
x = Yli^i Xib~ l and y = Yli^i Vib~ l , © is defined as 

oo 

x © y = z i b ~^ 

1=1 

where Zi = x» + yi (mod b). In the same manner, the operator denotes 
the digitwise subtraction modulo b. Further, we define a digitwise addition 
and digitwise subtraction for non-negative integers based on those base b 
representations. In case of vectors in [0, l) s or Nq, the operators © and 
are carried out componentwise. 

2.1 Walsh functions 

Walsh functions were first introduced in [33] for the case of base 2 and were 
generalized later, see for instance [5]. We first give the definition for the 
one-dimensional case. 

Definition 1 Let b > 2 be an integer and uj b = e 2m / h . We represent k G No 
in base b, k = kq + K\b + ■ ■ ■ + ft a _i& a_1 with Ki G {0, . . . , b — 1}. Then, the 
k-th b-adic Walsh function &walfc : [0, 1) — > {1, o;&, . . . , c^ -1 } is defined as 

for x G [0, 1) with base b representation x = Xib" 1 + x^" 2 + • • • , unique in 
the sense that infinitely many of the Xi are different from b — 1. 

This definition can be generalized to higher dimensions. 

Definition 2 For dimension s > 2, let ) G [0, l) s and k = 

(/ci, . . . , k 8 ) G Ng. We define 6 wal fc : [0, l) s -»• {1, u b , . . . , u^ 1 } by 

s 

b wal k (x) = JJ 6 wal^(xj). 

3=1 

Since we will always use Walsh functions in a fixed base b in the rest of this 
paper, we omit the subscript and simply write wal^ or wal^. 
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2.2 Polynomial lattice rules 

We introduce some notation first. For a prime b, let F& be the finite field 
containing b elements {0, . . . , b — 1} and by ¥b((x~ 1 )) we denote the field of 
formal Laurent series over F&. Every element of ¥b((x~ 1 )) is of the form 

L = ^tix~\ 

l=w 

where w is an arbitrary integer and all ti G F&. Further, we denote by F&[x] 
the set of all polynomials over Fj. For a given integer m, we define the map 
f m from ¥ b ((x~ 1 )) to the interval [0, 1) by 

5>-) = E w-- 

We often identify k G No, whose 6-adic expansion is given by k = k,q + nib + 
■ ■ - + K a _i6 a_1 , with the polynomial over ¥b[x] given by k(x) = kq + kix + - ■ ■ + 
K a _ix a_1 . For = (fci, . . . , k s ) G (F b [x]) s and q = (q u ...,q s ) G (F b [x]) s , we 
define the 'inner product' as 

s 

k q = ^kjqj G ¥ b [x], 
3=1 

and we write q = (mod p) if p divides q in Ff,[x]. 

The definition of a polynomial lattice rule is given as follows. 

Definition 3 Let b be prime and m, s G N. Let p G ¥b[x] be an irreducible 
polynomial with deg(p) = m and let q = (qi,...,q s ) G (F5[x]) s . Now we 
construct a point set consisting ofb m points in [0, l) s in the following way: For 
< h < b m , identify each h with a polynomial h(x) G Fj[z] o/deg(/i(x)) < m. 
Then the h-th point is obtained by setting 

The point set {xq, xi, . . . , Xb^-i} is called a polynomial lattice point set 
and a QMC rule using this point set is called a polynomial lattice rule with 
generating vector q and modulus p. 
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We add one more notation and introduce the concept of the so-called dual 
polynomial lattice of a polynomial lattice point set. For k £ N with 6-adic 
expansion k = Ko + Kib+- ■ ■ J r n a ^ib a ~ 1 , let tr m (/c) be the polynomial of degree 
at most m obtained by truncating the associated polynomial k(x) £ ¥b[x] as 



tr m (k) = k + Kix H h « m _ix m \ 

where we set K a = • • ■ = ft m _i = if a < m. For a vector k — (k\, . . . , k s ) £ 
Nq, we define tr m (fc) = (tr m (fci), . . . , tr m (k s )). With this notation, we in- 
troduce the following definition of the dual polynomial lattice D ± . 

Definition 4 The dual polynomial lattice of a polynomial lattice rule with 
modulus p £ ¥ b [x], deg(p) = m, and generating vector q £ (F 6 [a;]) s is given 
by 

D L = {keW : tr m (Jfe) • q = (mod p)}. 

The following important lemma relates the dual polynomial lattice to nu- 
merical integration of Walsh functions. 

Lemma 1 Let {xq, . . . } be a point set of polynomial lattice rule with 

modulus p £ Ffe[x] ; deg(p) = m, and generating vector q £ (F^fx]) 5 and let 
D be its dual polynomial lattice. Then we have 

b m — 1 

1 v — v , , s f 1 if k £ D- 



1 f 

— £ wal fc (^) = j 

i=0 ^ 



otherwise. 



Proof. This follows immediately from Definition [4j [E2 Lemma 10.6] and [T3| 
Lemma 4.75]. □ 

2.3 Owen's scrambling 

We now introduce Owen's scrambling algorithm. This procedure is best 
explained by using only one point x. We denote the point obtained after 
scrambling x by y. For ) £ [0, l) s , we denote the 6-adic 

expansion by 



x j,l . x j,2 . 
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for 1 < j < s. Let y = (yi, . . . ,y s ) G [0, l) s be the scrambled point whose 
6-adic expansion is represented by 

for 1 < j < s. Each coordinate t/,- is obtained by applying permutations 
to each digit of Xj. Here the permutation applied to Xj^ depends on Xjj 
for 1 < I < k - 1. In particular, ^ = ^-(a^i), y jj2 = TTj jXjtl {x jj2 ) , y ii3 = 
^■,^,1,^,2(^,3), and in general 

= 7l j,x^ 1 ,...,x jtk _ 1 (%j,k), 

where iTj tXj 1) ... ia; . fc _ 1 is a random permutation of {0, ... ,6 — 1}. We choose 
permutations with different indices mutually independent from each other 
where each permutation is chosen uniformly distributed. Then, as shown in 
[271 Proposition 2], the scrambled point y is uniformly distributed in [0, l) s . 

In order to simplify the notation, we denote by IT,- the set of permutations 
associated to Xj, that is, 

n i = {^j,x j ^...,x jtk _ 1 ■ k e N, x jt i, . . . , Xj, k -i e {0, . . . , b - 1}}, 

and let II = (n 1; . . . , U s ). We simply write y = Tl(x) when y is obtained by 
applying Owen's scrambling to x using the permutations in II. 

2.4 Higher order digital nets 

Quasi-Monte Carlo rules based on higher order digital nets exploit the smooth- 
ness of the integrand so that they can achieve the optimal order of conver- 
gence of the integration error for functions with smoothness a G N. The 
result is based on a bound on the decay of the Walsh coefficients of smooth 
functions [7J. We refer readers to [8] for a brief introduction of the central 
ideas. Explicit constructions of higher order digital nets and sequences were 
given in [7J. 

There is also a component-by-component construction algorithm of higher 
order polynomial lattice rules [3]. Higher order polynomial lattice rules can 
be obtained in the following way. In Definition [3], we set p with deg(p) = 
n > m and replace v m with v n for the mapping function. Then a higher 
order polynomial lattice point set consists of the first b m points of a classical 
polynomial lattice point set with b n points (where n = am for integrands 
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of smoothness a). The existence of higher order polynomial lattice rules 
achieving the optimal order of convergence was established in [TJ] and the 
CBC construction was proven to achieve the optimal order of convergence in 
[3] . However, we have no generalization of the scrambling algorithm to higher 
order polynomial lattice rules which preserves the higher order structure. To 
work around this problem we use a different approach in this paper. Namely, 
we use the approach from 0, E] based on the interlacing of digital nets or 
sequences. 

We describe the interlacing algorithm in more detail in the following. 
Since the interlacing is applied to each point separately, we use just one 
point to describe the procedure. Let z G [0, l) ds , with z = (z±, . . . , z ds ) and 
consider the 6-adic expansion of each coordinate 

unique in the sense that infinitely many digits are different from 6 — 1. We 
obtain a point x G [0, l) s by interlacing the digits of d components of z in 
the following way: Let ), where 

oo d 
a=l r=l 

for 1 < j < s. We denote this mapping by D d : [0, l) d — > [0, 1) and we simply 
write Xj = T>d{z(j-\)d+i, ■ ■ ■ , Zjd)- Further we write 

x = V d {z) := (£>0i, . . .,z d ),V(z d+1 , . . .,z 2d ), . . . ,V(z {s _ 1)d+1 , . . . , z sd )) 

when x is obtained by interlacing the components of z. 

Order-<i scrambling for z G [0, l) ds proceeds as follows. Let II = (III, . . . , H ds ) 
be a random set of permutations. Then, the order-rf scrambled point y G 
[0, l) s is given by 

y = V d (n(z)). 

Hence, as stated in the previous section, we first apply Owen's scrambling 
to z G [0, 1)^ and then interlace the digits of the resulting point to obtain 
the point y. Again, we choose permutations with different indices mutually 
independent from each other where each permutation is chosen with the same 
probability. Then, as shown in [HI Proposition 5], the order-c? scrambled point 
y is uniformly distributed in [0, l) s . 

In this paper, we are interested in the use of polynomial lattice rules to 
generate a point set in [0, l) ds . For clarity, we give the definition of interlaced 
scrambled polynomial lattice rules below. 
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Definition 5 Let b be prime and m,s,d G N. Let p G F^x] be an irreducible 
polynomial with deg(p) = m and let q = (qi, . . . ,qa s ) G (F fe [x]) ds . Now we 
construct a point set consisting of the b m points in [0, l) s . For < h < b m , 
the h-th point is obtained by setting 

_( f h(x) gi (x) \ f h(x) q ds (x) \\ ds 

Zh -{ V ™{ p{x) - p{x) JJe[0,l) . 

Then let 

y h = V d (U(z h )), 

where the permutations II are chosen independently and uniformly distributed. 
We call {y , y 1 , . . . , y bm _i\ an interlaced scrambled polynomial lattice point 
set (of order d) and a QMC rule using the point set {y , y l: . . . , y^^} an 
interlaced scrambled polynomial lattice rule (of order d). 



3 Variance of the estimator 

We consider the following Walsh series expansion for / G -^([O, l] s ) 

f( x ) ~ Yl f( k )™ al k{x), 

fceNg 

where the Walsh coefficients f{k) are given by 

f(k) = / f(x)wal h (x)dx. 

In the following, let {y , . . . , y^^} where y { G [0, l) s be an interlaced 
scrambled polynomial lattice point set. Further, we denote by {z , . . . , z^n^} 
where Zj G [0, l) ds the polynomial lattice point set with ds components with 
modulus p G Fb[x] an irreducible polynomial with deg(p) = m and with 
generating vector q = (q 1 , . . . , q^) G (F 6 [x]) d,s . We approximate the integral 

J (/) = f[a,i]° f( x ) dx b y 

6 m -l 

hf) = w E fly*)- 

i=0 
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Since the order-rf scrambled point is uniformly distributed in [0, l) s as shown 
in P Proposition 5], this estimator is unbiased, that is, E[J(/)] = /(/). It 
follows that the mean square error equals the variance of the estimator. Thus, 
in the following, we concentrate on the variance of the estimator denoted by 



/(/) - E[7(/)] ■ 



V[I{f)] = E 

The following notation is needed for deriving the lemma below. Let I = 
(h, . . . , l s ) G N^ where lj = (l {j _ 1)d+1 , . . . , l jd ) and 

B d , hs = {(h, k ds ) G N d s : L^-'J < % < b l * for 1 < j < ds}. 

In an analogous manner to T> d , we define a digital interlacing function £ d for 
non-negative integers. For k±, . . . , k ds G No, we represent the 6-adic expansion 
of kj by kj = Kj j0 + Kjflb + • • • for 1 < j < ds, where Kj tU G {0, 1, . . . , b — 1} 
(and where Kj >u = for all u large enough). Then, £ d denotes the following 
mapping from [k\, . . . , k ds ) G Nq s to (k[, . . . ,k' s ) G Nq, where 



iy—1+ad 

a=0 r=l 



for 1 < j < s. Then we define the following sum of the Walsh coefficients of 
/ over k G B d j, s , 

°iaf)= E i/(^(*))i a . 

and we introduce 

^ fe m -l (is 

rUg,p) = ^ E J]e [wau,(n,( 2ij ) e n.fe))] , 

i,i'=0 j=l 

where k = (hi, k 2 , . . . , k ds ) G B^i^ is an arbitrary element. We note that 
Ti td (q,p) is independent of the choice of k G i3d,« jS , and depends only on the 
point set {z , . . . , z&m_i}, see [9]. According to [HI Lemma 7], we have 

v(i(f))= £ < M (/)r M (g,p). (2) 

JeNjf \{o} 

By applying the property of polynomial lattice rules to this expression of 
V(I(f)), we obtain the following lemma. 
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Lemma 2 Let d > 1 and f G L 2 ([0, 1} S ). Let the estimator I be given by 

j 6 m -l 

A/) = 

i=0 

where {y , . . . , y&m^} an interlaced scrambled polynomial lattice point set 
with generating vector q and modulus p. Then, we have 

mu)) - E E E i. P) 

where Ids — {1; • • • , ds} ; |Z u |i = ^ Jgn ^- and -D -1 iae dual polynomial lattice 
for the polynomial lattice point set with generating vector q and modulus p. 

Proof. This follows immediately from [T3| Corollary 13.7]. □ 



3.1 A bound on the Walsh coefficients 

Below we define a variation V^(/) of order a > 1 for functions / : [0, l] s — > 
R. See [HI pp. 1386, 1387] for a derivation of this definition. In particular, 
in [9] it is shown that if the partial derivatives g^i dx U are continuous for 
a given . . . , a s ) G {1, . . . , a} s , then 

vi s \f) = U 

To define the variation Va , let J = Y\'*t 1 [(ijb~ l: > , (aj + l)b~ lj ), where 
< aj < b l i and lj G N for 1 < j < as. The set V a (J) = {V(x) : x G J} is 
the product of a union of intervals except for a countable number of points 
(see [§]). Let a G {1, ... , For t G [0, l) s and cci, . . . , x s G (-1, l) s we 
define the difference operator 

A a (t; Xl ,...,x s )f= Yl E (-1) M+ - +M 

DiC{l,...,ai} v s C{l,...,a a } 

x / I h + E •■•>*« + E 



0oi 
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da; 
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Then we define the generalized Vitali variation 



V°(f) 




Vol(£> Q (J))sup 



A a (t; xx,.. -,x s )f 



Tlj=l riril X j,r 



1/2 



where the first supremum sup-p is over all partitions of [0, l) as into subcubes 
of the form J = U%ii a jb~ lj » ( a j + l)b~ lj ) with < aj < & and lj £ N for 
1 < j < Ois, and the second supremum is taken over all t £ T> a (J) and cc 3 - = 
• • • , ^j,aj) with Xj >r = Tj tr b~ a ^ j ~ 1 ^~' r where r^ r £ {1 — b, . . . , b — 1} \ {0} 
for 1 < r < ttj and 1 < j < s and such that for all the points at which / is 

evaluated in A a (t; x u . . . , x 8 ) are in I> a (n^i[ 6 ~ ii+1 L a i/ fo J > & ~' J+1 K'/ & J 

For ^ it C {1, . . . , s} let \u\ denote the number of elements in u and 

let Va (f u ',u) denote the generalized Vitali variation with coefficient a £ 
{1, . . . , a}'"' of the | w | -dimensional function 



[0,1] S -I'"I 



f(x) das{i i ... ia }\ u . 



For u = we set f% = f, Q1 u f(x)dx and we define Va\h',^) = \f$\- Now 
we define the generalized weighted Hardy and Krause variation of / of order 
a by (cf. p. 1387]) 



v a (f)=[ E ^ E (^ |u|) (/u;«)) sl 

y«C{l,...,s} ae{l,...,a}l»l 

where (■y u )ueu is a sequence of nonnegative real numbers and U = {u C N : 
|w| < oo}. 

Let / : [0, l] s ->■ R and let 

f( x ) = E 9u(x u ) 

uC{l,...,s} 

denote the ANOVA decomposition of /, that is, g% = Jj ^ s f(x) dx and 
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where v C u means that v is a proper subset of u. We have g u (x u ) dxj = 
for j G u and |^ = for j (£u. Then we have 

W=7u 1/ \ M (W«). 

Let / = (h, . . . , G Nq s and let u = {i G {1, . . . , ds} : lj > 0}. Then we 
denote I by (Z U) 0). Let v(u) C {1, . . . , s} denote the set of 1 < i < s such 
that u n {(i - l)d + 1, (i - l)d + 2, . . . , id} ^ 0. Then 

Vd,(i u ,o),s(9v) = if w ^ u(m). 

Thus we have 

&d,(lu,0),s(f) = a d,(lu,0),s(9v(u))- 

Using [9], Lemma 9], we therefore obtain 

* d ,(i u ,o)M) = ^(i u , ),s)(9v(u)) < 2 l ^ )|max(d - a ' 0) /3(/ u ,0) v ^K,(/), 

where the definition of /3(l v ,0) is given as follows. Let itj = u n {(z — 
+ 1, and a, = min(a, |itj|) for 1 < z < s. Let /3j = (6 — 

for j G Ui and 1 < i < s. Let p i>x < /3 ij2 < ■ • • < 
for 1 < z < s be such that {/^i, . . . , A.a,} = {/^j : j S u^}, that is 
{Aj : 1 < J < «i} is just a reordering of the elements of the set : j G itj}. 
We define /3(l u ,0) as 

/3(^o)=nn/ 3 M- 

i=i j=i 

In the following lemma we provide a bound on the coefficients (3(l u , 0). 

Lemma 3 Let a, d, s G N. For / = (Z 1} . . . , £ ds ) G N$ s Zet /3(Z U , 0) be given as 
above. Further, we define fi(k) = a — 1 fork G N whose b-adic representation 
is given by k = k,q + + • • • + n a -ib a ~ l where n a -i 7^ 0. T/ien we /ia^e 

0) < Y[b~ min{a ' dMk ^. 

for all k G B d ^ ufi ), s - 
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Proof. First, we consider the case with d < a. Since we always have |itj| < 
d < a in this case, ccj = |itj|. Thus, from the definition of 0), we have 

s 

0(iu,o) = Y[Y[{b - i)b~ j+ ^ i)d ~ {i ^ i)d 

i=l j(zUi 

<{b- i)i u ir |u| J]V M( ^ )d 



< 



jeu 



Next, we consider the case where d > a. We denote by l'^^ d+ i, ■ ■ ■ ,l' id & 
descending reordering of l(i-i)d+i, ■ ■ ■ , kd, that is, Z[ i _ 1 - )d+1 > ■ • • > l' id where 
{i[i-i) d+ i, l'id} = 0(i-i)*H> • • • , kd}- Then we have 



P(iu,o)<H(b-irb- a >j[b-u 

i=l i=l 



(j-ljd+j 1 J 



< JJft-^^iKViw), (5) 



i=l 



where ft'c*- 1 )^ 1 < < b 1 ^- 1 ^ for 1 < z < s and 1 < j < d. If 

\ui\ < a < d, then «j = |itj| and it is obvious that we have 

1 a% 1 

for 1 < z < s. Otherwise if a < |itj| < d, then ai = a and, given that 
J(i-i)d + i> • • • > l'id is a descending reordering of • • • , ^d, we have 

j=l l| jr' = l je« 4 

for 1 < i < s. Thus, we have for both of the cases 

j=i jeu % 
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By inserting this inequality into (j5J), we have 

s 

/0(l«,O) < JJft- '^^ 

i=l 

The result follows by combining this inequality and (J4]). □ 
3.2 A bound on the variance 

Using Lemma [2] and the bound on the Walsh coefficients given in the previous 
subsection, we can now prove a bound on the variance of the estimator. We 
can then use this bound to introduce a quality criterion for the construction 
of interlaced scrambled polynomial lattice rules. 
Let us define 

(lr\ - I 1 If * = 0, 

r <*A K ) - S j-(2mta(a,d)+i)M(fe) otherwise. 

For k = (h,...,k ds ) E Nq s , let r(fe) = Uf^rfc) and r(k u ) = U jeu r ( k j) 
for k u = (kj)j £u . The following corollary gives a bound on the variance of 
the estimator. 

Corollary 1 Let a,d E N. lei / : [0, l] s -> R safe/y 1/ Q (/) < oo. Let the 
estimator I be given by 

n=0 

where {y , . . . , y^m^} is an interlaced scrambled polynomial lattice point set 
of order d > 1 w't/i generating vector q e (Ffe[x]) ds and modulus p. Then we 
have 

V(I(f))<V*(f) 41^)1^-^7^) £ r Q , d (fc u ,0). 

(fc^.ojei?- 1 - 

where D 1 - is a dual polynomial lattice of the polynomial lattice rule with 
generating vector q and modulus p as in Definition^ and v(u) C {1, . . . , s} 
is the set of alii G {1, . . . , s} such that un{(i—l)d+l, (i—l)d+2, . . . , id} 7^ 0. 
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Proof. From the bound on the Walsh coefficients given in the previous sub- 
section and Lemma [31 we have 



°i,o) A .(/) < 4 M " )|max(d - a ' 0) 7„ W ^ 2 (/)n^ 



-2 min(a, d)fi(kj ) 

for all (k u , 0) G Bd,(i U) o),a- Inserting this inequality into (EJ, we have 



a d,{i u ,o),s(f) 



^(A/)) E E E (6_ i)H6li»|i-H 



(fc u ,0)GD ± 

(fc„,0)GD- L 

Hence, the result follows. □ 

Remark 1 ,4s can fre seen from the proof, even if we replace the definition 
of r a ^(k) for k > with 

we s£i/£ reacn tae same form with Corollary Q| giving a tighter bound on the 
variance except the case with base 6 = 2. By dropping (b — 1) _1 from the 
above expression, however, we can obtain the result also for d > a, as shown 
later in Remark [H This implies the robustness of the constructed rules in 
terms of the smoothness of function. 

We denote the double sum in Corollary [T] by 

B aAj (q,p)= A^™< d -^ Mu) E W(*^0). (6) 

This value depends on the smoothness a, the weights (j u )uci 3 , both of which 
come from the function space, the interlacing factor d and the polynomial 
lattice rule with ds components. We note that it is independent of a partic- 
ular function /. Thus, it is possible to use B aAj (q,p) as a quality criterion 
for the construction of interlaced scrambled polynomial lattice rules. The 
following lemma gives a more computable form of B a ^y{q,p). 
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Lemma 4 For z E [0, 1), let 

b _ 1 _ ^2min(a,d)Llog i) 2j ^2min(a,d)+l _ ^ 



0a,dO) 



J J-2min(a,i 



where we set fr 2min (". rf )U°g,,oj _ q i e f B a>dn {q,p) be given by Then, we 
have 



b m -l 



1 + \\ (l + 0a,d(Zi,(j-l)d+fc)) 



fe=i 



i=o o^vci s jew 
In particular, for product weights j v = Yljev^j' we have 

1 b m -l s 

B a , dn (q,p) = -l + ^EIl[ 1 - 4 max(d - Q ' 0) 7, 

i=o i=i 
d 

fc=i 

Proof. Using Lemma [T], we can rewrite §5§ as follows 

1 6 m -l 



6" 

i=0 0^uC/ 4 fc„eNl»l 



6 m -l 



^E E 4l"""l-C--»'7„ w n 

i=0 0^«C/ ds jgu 



fc, = l 



(7) 

Following along similar lines as in the proof of [U Theorem 7.3], we have 

oo oo 1 b' — 1 

fc=l i=l A-^foi- 1 



^2 min(a,d)+ 



^ J _ ^2min(a,d)Llog b zJ ^2min(«,d)+l ]^ 

fr(fr2min(a,d) _ fj 
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We now rearrange ([7]). For a given v C J s we consider now sets w C I ds such 
that w(t>) = v. Thus we obtain 

d 

"I + JJ (l + 0a,d(>i,O-l)d+fc)) 
fe=l 

In case of the product weights 7^ = Yljev I31 t ne ^ as ^ expression can be further 
simplified to 

i=o o^vci s jeu 

i=0 j=l 
d 

+ 4 max ( d - a '°) 7 , n (1 + <p a , d (zi,u-i)d + k)) . 
k=i 

Hence the result follows. □ 



d 

-1 + 11(1 + ^(^0) 

fe=l 



4 Component-by-component construction of 
polynomial lattice rules 

4.1 Convergence rate of the variance 

The CBC construction algorithm was first introduced by Korobov [16], and 
reinvented later by Sloan and Reztsov [30J to construct a generating vector of 
lattice rules. The same approach can be applied to polynomial lattice rules. 
In the following, we restrict qj, 1 < j < ds, to non-zero polynomials over F& 
with its degree less than m. Without loss of generality we set q x = 1. We 
denote by Rb, m the set of all non-zero polynomials over with degree less 
than m, i.e., 

Rb,m = ¥ b [x] : deg(g) < m and q ^ 0}. 

We note that \Rb, m \ = b m — 1. Further, we write q ro = (qi, ■ ■ ■ ,q ro ) f° r 
1 < r o < ds and Ij = {1, . . . , jo} for 1 < j < s. The CBC construction 
intended for this study proceeds as follows. 
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Algorithm 1 For a prime base b, a dimension s, an interlacing factor d, 
and integer m > 1 and weights 7 = ('y u )uci s ■' 

1. Choose an irreducible polynomial p G ¥b[x] with deg(p) = m. 

2. Setql = 1. 

3. For r = 2,...,ds, find q* by minimizing B Qt d t -y((q* _i,q ro ),p) as a 
function of q TQ G Rb, m ■ 

In the proof of Theorem [T] and its subsequent remark, we will use Jensen's 
inequality, which states that for a sequence (a^) of non-negative real numbers 
we have 

A 



for any < A < 1. 



Theorem 1 Let b be a prime and p G F&[x] be irreducible with deg(p) = m. 
Suppose that q* = (<?*, . . . ,q^ s ) is constructed using Algorithm^ Then, for 
all ro = 1, ... ,ds we have 



B a4a {q*p)<{b m -iy- 



O^«CJ 3 - _i uC7j _! 



/or all 1/(2 min(a, d) + 1) < A < 1, where r = (jo — l)d + do (jo, d G N and 
< d < dj, 

D a , A = -l + (l + C w ,A) a , 

and 

A 



ChaAX = max 



^_max(d— «,0) l)\ ^Amax((i-a,0)^ -p 



^ _ -2min(a,d) / ' 1 _ W-(2 min(a,d)+l)A 



Proof. We prove the result by following along the same lines as in the proof 
of [TQl Theorem 4.4]. We proceed by induction. First for r = 1, that is, for 
jo = 1 and d = 1, we have 



B aAl (l,p) = 4 m -( d - Q .°) 7{1} r a ,d{k) 



k=l 
b m \k 
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(2 min(a,d)+l)^(fcj) 



k=l 

b m \k 



4 max(d-a,0) 7{i} _ ^ b l-m b - 



(2min(o,d)+l)Z 



l=m 



^max(d— a,0) 



6-1 



7{i} 

^(2 min(a,d)+l)m J ^-2 mm(ct,i 



< (6 m - 1) 



^max(d— a,0)A^A 



b-1 



^ ^— 2m'm(a,d) 

for all 1/(2 min(a, d) + 1) < A < 1. Consequently, we obtain 

iWM<(fr m -i)^ Hi}Di,xy. 

Next, suppose that for some r = (jo — l)d+do (jo, do EN and < do < c?) 
we have qr* e i?r° which satisfies 



iWg;,p)<(& m -i) 



"I 



We denote tq + 1 = (ji — i)d + d\ (ji,d\ € N and < d\ < d). It is obvious 
that we have 



(ji,di 



(jo + 1, 1) if d = d, 
(jO) do + 1) otherwise. 



Now we consider 



(fe. u ,0)GD J 



) I max((i— a,0) 



(fc tt ,o)eZ) J 



uGI r 



7v(uU{r +l}) E r (^ ; ^ro+l)0) 

(fc Il ,fc ro+ i)GNl u l+ 1 
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=B aAl (q* ro ,P) +%r +l), 

where we define 

6{q T0+l ) := 4 KuU{rt,+1})|max( ' i - cl ' 0) 7,(«u { ,-o + i}) E r ( fc «> ^o+i, 0). 
«cir Q (fc tl ,fc ro+ i)eNl' 1 l+ 1 

(fc I1 ,fc ro+1 ,o)eD- L 

In order to minimize B a>d> ^((q* o , q ro +i), p) as a function q ro+ i, we only need 
to consider 9(q ro +i). Based on the averaging argument, that is, the minimum 
value of 9(q rQ+ i) is less than or equal to the average value of 9(q ro+ i), we have 
for 1 /(2 min(a, d) + 1) < A < 1 



6- 



K + i)<^r E ^ 



ro+1; 

< _ A A|«(MU{r-o+l})|max(£i-a,0) /v A 
-ftm_J 7^(«U{r +l}) 

(fc u ,fcr 0+ i)eNl"l+ 1 
(fc u ,fc ro+ i,o)eD ± 

^ A|u(nU{ro+l})| max(d— a,0)-,A 

E* MmU{t-q + 1}) 
6 m - 1 

MC7 rQ 

E E r\k u ,k ro+1 ,0), (8) 

q ro + l&R b , m (fc u ,fc ro + 1 )GNl' 1 l + 1 
(fc u ,fc ro+1 ,0)GD ± 

where we have used Jensen's inequality in the second inequality. For a fixed 
u Q Ir of the outermost sum in (EI), if K +i is a multiple of b m , we always 
have tr m (/c ro+ i) = and the corresponding term becomes independent of 
q ro +i, or otherwise we have tr m (Av +i) ^ and tr m (/c ro+ i)g ro+ i cannot be a 
multiple of p by considering that p is irreducible. Hence we have 

j^rzrj E E Ak u ,k ro+1 ,o) 

<?r + l€Ri,,m (fe„,fc ro + 1 )GNl"l + 1 

(fc li ,fc ro+ i,o)e-D ± 

oo 

= E rA (^+i) E rX ^ 

b m \k r()+1 tr m (fe„)-q„=0 (mod p) 
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1 oo 

fer +i=i fc„eNl u l 

6 m t*V + i tr m (fc„)-g u ^0 (mod p) 

For the first term of the right-hand side in (Q, we have 



E ^(k n+1 ) = E 6 ~ 



(2 min(Q,c/)+l)A/i(fc ro _|_i) 



+ 1 — 1 + 1 — 1 

b m \k rQ + 1 b m \kr +l 

b-1 



E* (1 - 



(2mm(a,d)+l)\)l 



b m 

l=m 



For the second term, on the other hand, we have 



ftm _ J E rA (^0+l) 

fcr + l = l 

1 m-1 oo 

—EE ^w+^jE E 



l ~° k r() + 1 =b l *— '"■ K ro + l 



m—1 

-(2min(a,d)+l)AZ 



L_^(&-i)^- 



6 m 

z=o 



(2min(a,d)+l)AZ 

6 m 



-^^(6- l)(6 m - l)6'- m 6 

Z=m 

_^ _ £(l-(2min(a,d)+l)A)Z , & £ ^(l-(2 min(a,d)+l)A)Z 

1=0 l=m 

By inserting these equalities into ([9]), we have 



~I E E r x (k u ,k ro+ i,0) 



b m 

9r + l6R 6 , m (fc li ,fc ro + i)GNl u l + 1 
(feu ,AVq + 1 i 



6-1 



y^ & (l-(2mln(a,<Q+l)A)i ^ r A (fc u ) 



6" 

2=m feuGNl 1 
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6-1 



m—1 



(2mra(a,d)+l)A)Z 



E 



r\k u ) 



1=0 



tr m (fc„) <j 11 ^0 (mod p) 



< _^ L 7 ) ( 1 ~( 2 min(a,d)+l)A)/ 



«=0 



E rA ( fc «) 

fc„eNH 



6-1 



(fcm _!)(]_ _ 61-(2mm(a,cO+l)A) 



n 



E^ 



Here the sum in the product is given by 



E rA (^) = E b ~ 



(2mm(a,d)+l)Xfj,(kj) 



;g^-(2min(a,d)+l)AZ 

Z=0 

6-1 

^ _ ^l-(2min(a,d)+l)A ' 



Thus, from (IH1) we obtain 



uCI rn 



(10) 



Recall that tq = (ji — l)d + d\ — 1. Let Ji := {1, . . . , (ji — l)c/} and J2 '■— 
{(ji — l)d + 1, . . . , (Ji — l)d + di — 1}. In case of di = 1, the set J 2 is set to 
be empty. Every subset u C I ro can be classified by whether it includes some 
element of J\ and by whether it includes some element of J2. Since {r$ + 1} 
is one of d components for the ji-th dimension, whether or not u includes 
some element of J2 does not affect J v (uu{r +i})- From this observation, we 
have 



Ea ^M+i _ H n2 l +1 \^ a r 

7v(«U{r +l}) o 6,a,(i,A — 2—i b,a,d,X 2—i 'v(«l)U{j'i} ' 
uCI rQ U2CJ2 U1CJ1 



ui\ 

b,a,d,X' 



By considering the terms associated with a certain u (u C in the inner 

sum, at least one component of {(j — l)c? + 1, . . . , j<i} for all j G w must be 
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chosen. Thus, 



7^(nU{r + l})°b,a,d,A 

U<ZI ro 



--C b ,a,d,\(l + C b ^ aA x) dl 1 E 7«u{ji} II ( _1 + (1 + ^b,a,d ) X) d ) 

z Ch,a,d,\(l + Cb^xY 1 ' 1 E 7uU{ji} ( _1 + (1 + Cfe,a,rf,A) d )' 

uC7 3 - 1 _ 1 

- (Dd l7 x — -D<Zi-i,a) J2 TiuOi}^S- 



Finally, we have 

5a,d, 7 « + i,:p) = B aAj (q* ro ,p) + 



<(6 m - 



+ (6 m - 1)- 



1 l/A 



(A^a - A^-i.a) E riu{h} Dl £ 



uCJ. 



In case of < do < d, we have ji = jo an d rfi = d + 1. Using Jensen' 
inequality, we obtain 



B a4n (q* ro+1 ,p)<(b m -l)-> 



E ^ D T + D ^ E tW#a 



In case of <io = d, we have ji = jo + 1 and d± — 1. Again by using Jensen' 
inequality, we obtain 



(9*o+i' P) 
<(6 m - 1)"* 



V t a z) |u [ +1 

/ „ lu u d,\ 



+ (6 m - 



"I 



uCI 



<(b m - 



E 7^i +1 + A^ E -rWt 1 



Jl-1 

1 

A 
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Hence, the proof is complete. 



□ 



Remark 2 We have shown that we can construct an interlaced polynomial 
lattice rule which satisfies 

B a , dn (q*,p) < A b , s , a4n , x b-^^ d ^ m+s , 

l + 2a' 

for all 5 > 0. For a < a < d and 7=7 1+2a , it follows from Jensen's 
inequality that 

{k u ,0)£D ± 

= E ^\uT E C/ a ( fc -°) 

(fc li ,0)GD- L 

/ * x l + 2a' 

< B aAl (q*,p) 

— b,s,a,d,~f,X 

for all 5 > 0. 27ms means that interlaced polynomial lattice rules constructed 
component-by- component for functions of smoothness a using an interlac- 
ing factor of d still achieve the optimal rate of convergence for functions of 
smoothness a' as long as a < a' < d holds. Our observation is similar to 
that of the classical polynomial lattice rule shown by fMj, while we note that 
it is opposite from propagation rules JB, Theorem 3.3} which states that a 
higher order net which achieves an optimal rate of convergence for function 
of smoothness a can achieve an optimal rate of convergence for function of 
smoothness a' for all 1 < a' < a. 

We briefly discuss the tractability properties of the interlaced scrambled 
polynomial lattice rules. 

Corollary 2 Let b be a prime base, p E ¥b[x] be irreducible with deg(p) = m. 
Suppose that q* is constructed according to Algorithm [H Then we have the 
following: 
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1. For general weights, assume that 



V <y x D ]u] 



< oo. 



0^uCI a 



Then B a ^ n {q* ,p) is bounded independently of the dimension for all 
l/(2mm(a,d) + 1) < A < 1. 

2. For general weights, assume that 



lim sup 



A n l u l 

d,X 



< OO. 



for some q > 0. Then the bound B a ^(q*,p) depends polynomially on 
the dimension with its degree at most q for all 1/(2 min(a, d) + 1) < 
A < 1. 

3. For product weights ^ u = Yljeulj' assume that 



< oo. 



Then, B ai d n {q* ,p) is bounded independently of the dimension for all 
l/(2min(ajd) + 1) < A < 1. 

4- For product weights ^ u = Ylj £u lj, assume that 



lim sup — 

s^oo log S 



< OO. 



for some q > 0. Then, the bound of B ol) d )1 {q* ,p) depends polynomially 
on the dimension for A = 1 . 

Proof. The results for general weights follow along similar lines as the proof 
of [HI Theorem 3]. In case of product weights 7„ = Ylj^Jj, we have from 
Theorem [1] 



B a>d Jq*,p)<{b m -l)~D 



A 



-i+ n(i+^D d>x ) 



The results for product weights follow along similar lines as the proof of [10| 
Corollary 4.5]. □ 



27 



4.2 Fast construction for product weights 

We now show how one can use the fast component-by-component construc- 
tion to find suitable polynomials qi, ■ ■ ■ ,qds £ ^b[x] of degree less than m for 
product weights. From Definition [5] and Lemma HI we have 

B a ,d,-f{{qi, ■ ■ -,qds),p) 

1 b m -l s 

= - 1 + f e n I 1 - ^-"h 

i=0 j=l 

+ 4^ 7J - n (i + ^ ))) ] • 

According to Algorithm (TJ set gi = 1 and construct the polynomials 
q2, ■ ■ ■ , qds inductively in the following way. Assume that q2, ■ ■ ■ , q ro are al- 
ready found. Let r = (jo - l)d+do andr +l = (ji — l)d+di (jo,do,ji,di G N 
and < d , d\ < d). As in the proof of Theorem [1], (ji, d\) = (j + 1, 1) if 
d = d, or otherwise (ji,dx) = (jo, do + 1). Here we introduce the following 
notation 



and 



for < i < 6 m . We note that <5i,. ro = 1 when di = 1 (or d = d). 

Since p is an irreducible polynomial over F&[a?], there exists a primitive 
polynomial g in ¥ b [x]/p, that is {g°(x) = g b " l ~ 1 (x) = 1, g 1 (x), . . . , g bm ~ 2 (x)} = 
(¥ b [x]/p) \ {0}. Using the above notation, we have 



B a ,d n ((q r o,9 z ),p) 

1 b m -l 

- ~ 1 + U-a E P ^o [l - 4 maX ^> ) 7jl 



i=0 
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+ ^-"'H^ro (l + <f>a4 (v m ( 9Z(X p 9 { 1 (X) 

for < z < b m , where g~ l (x) = g bm ~ 2 (x) is the multiplicative inverse 
of g(x) in ¥f,[x\/p (which is also primitive). The aim here is to compute 
B a ,d n ((q ro ,g z ),p) for < z < b m and choose z such that B a>dn ((q ro , g z °),p) < 
Ba,d,-y((Q ro i 9 Z ) iP) f° r all < 2 < b m . Since we only need to compare the 
values of B a dn ((q ro , g z ),p) for different values of z, we only need to compute 
the terms which depend on z. That is, it is sufficient to compute 



i=l 

We define the circulant matrix 



g z 

p(x) 



A = [<f>a,d ( V 

and a = (a±, . . . , a;,m_i) T with 



fit 



g z \x) 

p(x) 



Kz,i<6" 



Let b = Aa with b = (b\, ... , bb^-i) 1 ■ Then zq is the integer 1 < zo < b m 
which satisfies b Zo < b z for 1 < z < b m . Therefore we set g ro +i = g zo . 

Since the matrix A is circulant, the matrix vector multiplication Aa can 
be done using the fast Fourier transform as shown in [25||26]. Thus we obtain 
a fast computation of the vector b. 

After finding q ro+ i, Pi >ro and Qi )T0 for < i < b m are updated as follows. 
If di = d, 

Pi, ro +i= P vro+1 [l-4 max ( rf -«'°) 7jl 

Qi,r +1 = 1- 

Otherwise if < d\ < d, 



Pi,ro+l Pi,roi 

Qi,r +1 — Qi,r a (l + 4>a,d \ V- 



q TQ+1 {x)i{x) 
p(x) 
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Then, we proceed to the next component. Unlike for classical polynomial 
lattice rules, such as [2J [26J , here we are required to store not one but two 
vectors (Pi, ro ) and (Qi, ro ) in memory. By this slight increase of memory, the 
fast CBC construction using the fast Fourier transform can be applicable. 
The construction cost is of order 0(dsmb m ) operations using 0(b m ) memory. 
This compares favorably with the construction of deterministic higher order 
polynomial lattice point sets in [I] where the construction cost was of order 
0(sN a a log N) operations using 0(N a ) memory. 

5 Numerical experiments 

Finally, we present some numerical results for the bound on the variance of 
the estimator in weighted function spaces with smoothness of a > 1. In our 
computation, the prime base b is always fixed at 2 and the product weights 
are considered. In Figure [IH3l we show the values of B a ^ n (q,p) from m = 4 
to m = 16 with different choices of a, d, 7 and s, which are obtained by using 
the fast CBC construction. As proven in Theorem [TJ the CBC construction 
of interlaced scrambled polynomial lattice rules achieves the variance of the 
estimator of order n -( 2mm ( a < d )+ 1 )+ s (5 > 0), where n = 2 m in this experiments. 

In Figure [H s — 1 and 71 = 4- max ( d - a '°) with various choices of (a,d). 
When a — 1, the decay rate is of order n~ 3 regardless of the interlacing factor 
d. The difference is due to the coefficient which grows as d increases. As a 
and d increase simultaneously, the convergence rate is of orders around n~ 5 
and n~ 7 for d = 2 and d = 3, respectively. We note that in case of a > d, 
B a ,d,-y(q,p) is independent of a and that the resulting values of B a ^,-y(q,p) 
are equal to those for the case of a = d. 

Figure [2] shows the results for s = 2 and 71 = 72 = 4~ max ( d - Q '°) with 
various choices of (a,d). The orders of the convergence rate are about n -2 ' 5 , 
n -3,5 , and n -4 ' 5 , for a = 1, 2, and 3, respectively. 

In Figure [3] we compare the results for two different product weights with 
a — d — 2 from s — 1 to s — 5. One is jj = 1, the other is jj = j~ 2 . The 
latter implies a decreasing importance of the successive coordinates. It is 
obvious that as s increases the difference between two results becomes larger. 
Meanwhile, we can observe the asymptotic improvement of the convergence 
rate from the results for s — 5. 
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Figure 1: Values of B a ^ n (q,p) for s — 1 and 71 = 4- max ( d - Q >°) with various 
choices of (a, d) = (1, 1), (1, 2), (1, 3), (2, 2), (2, 3), (3, 3) marked respectively 
by red square, red circle, red triangle, blue down-pointing triangle, blue dia- 
mond, and green pentagon. 
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Figure 2: Values of B at d n (q,p) for s = 2 and 71 = 72 = 4 max ( d Q >°) with 
various choices of (a, d). The colors and marks used are the same as Figure 

m 
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Figure 3: Values of B a dn (q,p) from s = 1 to s = 5 with (a, d) = (2,2) for 
two different product weights: jj = 1 by blue and jj = j~ 2 by red. 
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